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THREE-DIMENSIONAL  THERMOELASTIC  PROBLEMS  OF  PLANES 


OF  DISCONTINUITIES  OR  CRACKS  IN  SOLIDS1 

by 

2  3 

M.  K.  Kassir  and  6.  C.  Si h 

Abstract 

Presented  in  this  paper  is  a  general  formulation  of  the  three- 
dimensional  thermoelastic  equations  (uncoupled)  for  problems  involving 
crack-like  imperfections  or  planes  of  discontinuities  of  some  kind  in 
solids.  Harmonic  functions  are  constructed  from  which  the  stresses  and 
displacements  may  be  obtained.  The  thermoelastic  potential  for  symmetric 
distribution  of  temperatures  on  the  crack  surfaces  is  related  to  Boussi- 
nesq's  three-dimensional  logarithmic  potential  for  a  disk  in  the  shape  of 
the  crack.  The  mass  density  of  the  disk  is  found  to  be  proportional  to 
the  prescribed  temperature  gradient  normal  to  the  crack  plane.  As  an 
application  of  the  theory,  closed  form  solutions,  in  terms  of  complete 
and  incomplete  elliptic  integrals  of  the  first  and  second  kind  are  given 
for  a  flat  elliptical  crack  whose  surfaces  are  exposed  to  uniform  temper¬ 
atures  and/or  temperature  gradients.  The  possibility  of  extending  the 
Griffith- Irwin  theory  of  fracture  to  cracks  in  thermal  environments  is 
also  discussed. 


^he  results  presented  in  this  paper  were  obtained  in  the  course  of  re¬ 
search  carried  out  under  Contract  Nonr  610(06)  with  the  Office  of  Naval 
Research,  United  States  Navy.  This  paper  will  be  published  in  the  Pro¬ 
ceedings  of  the  Third  Southeastern  Conference  on  Theoretical  and  Applied 
Mechanics. 

2 

Instructor  of  Mechanics,  Lehigh  University,  Bethlehem,  Pennsylvania. 

3 

Professor  of  Mechanics,  Lehigh  University,  Bethlehem,  Pennsylvania. 


Notation 
a,  b 
f,  9,  h 
k,  k' 

kj(j=l  ,2,3) 
q(x,y) 
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s,  t 
u,  v,  w 

VVwo 

X,  y,  z 
A,  B,  C 
E(k),  K(k) 

F,  G,  H 
I 

Q(s) 

Qo 

R,  Rq 
T(x,y,z) 

To 

Z 

a 

8 

Z,  n,  C 
x(s) 


Semi -axes  of  elliptical  crack. 

Harmonic  functions  of  x,  y,  z. 

Arguments  of  complete  elliptic  integrals, 
k  ^  k  —  1 . 

Stress-intensity  factors. 

Density  distribution. 

Constants . 

Distance  normal  to  crack  border. 

Dummy  variables. 

Rectangular  components  of  displacement. 

Displacements  on  crack  surface. 

Rectangular  coordinates. 

Multiplying  constants 

Complete  elliptic  integrals  of  second  and 
first  kind  associated  with  k,  respectively. 

Potential  functions  of  x,  y,  z. 

Integral . 

Function  of  s. 

Uniform  temperature  gradient. 

Length. 

Temperature  distribution. 

Constant  temperature. 


Coefficient  of  thermal  expansion. 
Angle  subtended  by  r  and  x-axis. 
Ellipsoidal  coordinates. 

Arbitrary  function. 
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y 

v 


p 


°xx»°yy»0zz 

Txy*Tyz*Tzx 

♦ 

♦ 

•m 

* 

|jl*  ,  |J/** 

A,  Q 


Shear  modulus. 

Poisson’s  ratio. 

Ratio  of  r  to  a. 

Normal  components  of  stress. 

Shear  components  of  stress. 

Angle  in  parametric  equation  of  ellipse. 
Vector  potential. 

Thermoelastic  potentials. 

Scalar  potentials. 

Functions  in  skew-symmetrical  problems. 


Introduction 

While  the  solution  of  problems  of  thermal  stress  for  thin  plates  and 
long  cylinders  has  received  considerable  attention  in  the  past4,  compara¬ 
tively  little  work  has  been  done  on  the  problems  of  three-dimensional 
theory  of  thermoelasticity  for  solids  containing  imperfections.  It  is 
known  that  when  a  temperature  field  is  disturbed  by  the  presence  of  cracks 
or  flaws,  there  is  high  elevation  of  the  local  temperature  gradient  accom¬ 
panied  by  thermal  stress.  Thermal  disturbances  of  this  kind  in  some 
cases  cause  crack  propagation  resulting  in  serious  damage  of  structural 
members.  Hence,  solutions  to  thermoelastic  problems  of  cracks  have  imme¬ 
diate  practical  value  in  fracture  mechanics. 


The  existing  literature  on  the  three-dimensional  aspects  of  thermal 
stress  around  cracks  has  been  limited  to  a  few  publications.  Moreover, 


4 

An  extensive  reference  of  papers  on  thermoelasticity  may  be  found  in 

[1,2]. 
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previous  works  are  confined  to  axially  symmetric  problems  dealinq  with  the 
"penny-shaped"  crack.  Such  a  restriction  is  mainly  due  to  the  mathematical 
difficulties  of  this  class  of  problems.  By  the  method  of  dual  integral 
equations  in  the  Hankel  transform,  Olesiak  and  Sneddon  [3]  investigated 
the  distribution  of  thermal  stress  in  the  neighborhood  of  a  penny-shaped 
crack  in  an  infinite  medium.  The  prescribed  temperatures  on  the  upper 
surface  of  the  crack  are  identical  with  those  on  the  lower  surface.  Using 
the  audition  of  symmetry  with  respect  to  the  crack  plane,  they  reduced 
the  problem  to  one  of  mixed  boundary  conditions  on  the  surface  of  a  semi- 
infinite  solid.  The  same  problem  was  also  treated  by  Deutsch  [4].  The 
case  when  the  thermal  conditions  are  applied  skew-symmetrical ly  to  the  crack 
surfaces  was  considered  by  Florence  and  Goodier  [5].  The  present  paper 
deals  with  a  class  of  thermoelastic  problems  with  cracks  of  more  general 
shapes  which  include  the  circular  cracks  in  [3-5]  as  special  cases.  It 
should  be  pointed  out  that  the  effects  of  both  inertia  and  coupling  between 
temperature  and  strain  fields  are  neglected  in  the  aforementioned  works. 

An  effective  way  of  solving  three-dimensional  boundary-value  problems 
is  to  construct  the  general  solution  in  terms  of  certain  arbitrary  func¬ 
tions  dictated,  in  part,  by  the  appropriate  field  equations  and,  in  part, 
by  the  topology  of  the  region  of  interest.  These  arbitrary  functions  are 
then  determined  such  that  the  boundary  conditions  of  the  problem  are  ful¬ 
filled.  This  is  precisely  the  approach  adopted  in  the  work  to  follow. 

With  the  aid  of  harmonic  functions,  a  general  solution  of  Navier's  equa¬ 
tion  including  thermal  effects  is  obtained.  By  restricting  the  analysis 
to  problems  with  cracks,  stress  and  displacement  expressions  are  derived 
for  two  types  of  problems,  one  of  which  concerns  with  stress  systems 
which  are  symmetrical  about  the  crack  plane  and  the  other  with  skew- 
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symmetrical  systems .  In  the  former  case,  the  steady-state  temperature  field 
for  an  arbitrary  region  of  exposure  (or  region  of  the  crack)  may  be  deter¬ 
mined  from  the  Newtonian  potential  for  a  disk  in  the  form  of  the  crack  plane. 
For  the  purpose  of  I1 .ustration,  detailed  solutions  in  ellipsoidal  coordi¬ 
nates  are  obtained  for  the  problem  of  heat  applied  to  the  surfaces  of  an 
“elliptically-shaped"  crack. 

Since  the  practical  aspect  of  this  paper  is  to  establish  a  criterion 
of  fracture  for  cracks  in  thermally  stressed  bodies,  an  examination  of 
the  stress  field  in  the  vicinity  of  the  crack  border  is  pertinent.  The 
local  stresses  are  found  to  have  the  same  functional  relationship  and 
inverse  square -root  singularity  as  those  obtained  in  isothermal  elastic 
bodies  subjected  to  surface  tractions  [6].  The  significance  of  this  result 
is  that  the  Griffith-Irwin  theory  of  fracture  [7],  oriqinally  developed  for 
bodies  maintained  at  constant  temperatures,  may  now  be  extended  to  predict 
the  onset  of  rapid  crack  propagation  caused  by  thermal  changes.  More 
specifically,  stress -intensity  factors  kj(j=l,2,3),  which  govern  the  sta¬ 
bility  behavior  of  cracks,  are  computed  and  shown  in  curves  for  the  cases 
of  thermal  conditions  applied  symmetrically  and  skew-symmetrical ly  to  the 
faces  of  an  elliptical  crack. 

Fundamental  Equations  of  Thermoelasticity 

When  the  influences  of  both  coupling  and  inertia  are  disregarded,  the 
general  thermal -stress  problem  separates  into  two  distinct  problems  to  be 
solved  consecutively.  The  first  is  a  problem  in  the  theory  of  steady- 
state  heat  conduction  which  requires  the  temperatures  T(x,y,z)  at  every 
point  of  the  body  to  satisfy  the  Laplace  equation  in  three-dimensions: 

v2T(x,y,z)  =  0  .  (1) 
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Once  the  temperature  distribution  has  been  found,  the  resulting  displace¬ 
ments  and  stresses  may  be  obtained,  respectively,  from  the  Navier's  dis¬ 
placement  equation  of  static  equilibrium 


1  2  „,l+v  ,  ^ 

JI2Z  77*“  +  7  “  =  2(l^a7T- 
and  the  Duhamel -Neumann  stress-displacement  relation 

a  =  y|vu  +  UV  +  [vV.U  -  (l+v)aT]  I J. 


(2) 


(3) 


in  which  u  is  the  displacement  vector,  o  is  the  stress  tensor  and  I  is  the 
isotropic  tensor.  The  gradient  and  Laplacian  operators  are  denoted  by  v 

p 

and  v  ,  while  y,  v,  and  a  designate  the  shear  modulus,  Poisson's  ratio, 
and  the  coefficient  of  linear  expansion  of  the  solid  whose  mechanical  and 
thermal  properties  are  assumed  to  be  isotropic  and  homogeneous. 

If  x,y,z  stand  for  the  Cartesian  coordinates,  the  solution  of  eq.  (2) 
for  three-dimensional  problems  with  geometric  discontinuities  on  the  plane 
z  =  0  takes  the  form 


U  -  <{>  +  ZVip*  ,  ip*  =  ip**  +  \p. 


(4) 


where 


v2<j> 


■  0,  72**  ■  0,  |jf  * 


and  the  vector  displacement  potential  $  has  components  <f>x ,  <j>y,  and  <j>z. 
In  eq.  (4),  ip**  is  the  scalar  displacement  potential  and  4/  the  thermo¬ 
elastic  potential.  Hence,  the  dilatation  becomes 


V.U  =  V .  <j)  +  “  . 


and  the  condition  of  equilibrium,  eq.  (2),  requires  that 


(5) 


(6) 
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(7) 


V.4  +  [  (3-4v)4>**  -  ip  ]  =  0. 

For  T(x,y,z)  =  0,  the  above  expressions  reduce  to  those  given  in  [6]. 

Using  eqs.  (3),  (4)  and  (6),  the  components  of  the  stress  tensor  o, 
in  terms  of  $  and  ip* ,  may  be  obtained: 


*  (♦+♦*>+  *  2z 

3X  Z  9z 


3X3Z 


(8b) 


'f-’k  <V*‘>+^+2z 


A* 

9y3Z  * 


(8c) 


For  the  purpose  of  discussing  symmetry  conditions  in  the  section  to  follow, 
the  other  stress  components  oxx,  oyy,  -rXy  are  not  of  immediate  interest. 


Thermoelastic  Problems  with  Planes  of  Discontinuities 

In  order  to  avoid  unnecessary  complication  in  the  analysis,  the  surfaces 
of  the  plane  of  discontinuity  or  crack  at  z  =  0*  are  taken  to  be  free  from 
applied  normal  and  shear  stresses,  since  the  method  of  solution  when  there 
are  such  applied  stresses  has  already  been  treated  in  [6,8].  Moreover,  the 
complete  stress  solution  can  be  split  uo  into  two  parts,  which  may  be  con¬ 
sidered  independently.  The  first  part  deals  with  the  application  of  surface 
temperatures  that  are  the  same  on  both  sides  of  the  crack  and  the  second 
considers  the  case  where  the  temperatures  on  the  upper  surface,  z  =  0+,  of 
the  crack  are  equal  and  opposite  to  those  on  the  lower  surface,  z  =  0  . 

The  corresponding  stress  systems  will  be  referred  to  as  the  symmetrical  part 
and  skew-symnetrical  part,  respectively. 


In  the  symmetrical  part  of  the  problem,  the  stresses  and  displacements, 
induced  by  thermal  changes,  will  depend  upon  the  variable  z  as  follows: 
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;  even  in  z 


(9a) 


u’  v’  axx’  °yy *  CTzz 


w,  iX2,  Ty; 


;  odd  in  z 


(9b) 


These  functions  are  required  to  be  continuous  outside  the  crack  region  on 

the  plane  z  =  0  and  the  odd  ones  shown  in  eq.  (9b)  must  be  zero  on  such  a 

plane.  Furthermore,  if  the  crack  surfaces  are  free  from  applied  stress, 

then  a  =  0.  In  view  of  symmetry,  the  condition  t  =  t  =0  must  hold 

xz  yz 

everywhere  on  the  plane  z  =  0.  Hence,  eqs.  (8b)  and  (8c)  give 


3 

3X 


(<t>z+<p*)  +  =  0 


z  =  0. 


(10a) 


(*  +^*)  +  =  0,  z  =  0. 

ay  z  dz 

At  this  point,  it  is  convenient  to  introduce  a  harmonic  funct4on  f(x,y,z) 
such  that  both  eqs.  (7)  and  (10)  are  satisfied,  i.e.,  by  letting 

*x  =  (1‘2v>  lx  +  fz  f£  dz  • 

4>  =  d-2v)  fy  +  /  dz  . 
y  z 

...  .  3f  3f 

♦z  =  -2(l-v)  aj  ,  ip**  =  ^  . 


(10b) 


(Ha) 

(Hb) 

(lie) 


where 

v2f(x,y,z)  =  0  . 


The  limits  of  the  integrals  in  eqs.  (11a)  and  (lib)  have  been  chosen  to 
satisfy  the  condition  of  regularity  of  the  displacement  at  infinity.  It 
follows  that  the  displacements  and  stresses  in  the  solid  may  be  expressed 
in  terms  of  two  real  harmonic  functions  f(x,y,z)  and  i|i(x,y,z).  They  are 
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u  *  n*  /“Hdz  + 

dz  + 

« =  -2d-v)  n  +  2  f  . 


3F 

3X  * 

(12a) 

3F 

m 

(12b) 

ay 

(12c) 

gxx  _  32f 
2 y  ~  ax7 


+  2v 


0  *2r 
J£  =  i4  + 
2y  ay7 


2v 


(13a) 


(13b) 


°zz  .  aF  ,  ,  a2F 

?r  ~ 9?  +  Z  15? 


(13c) 


lii.  .  (l-2u)  »!f_  +  f”  8?1_  dz  +  z  £L. 
2y  3x3y  z  3x3y  3x3y 


(13d) 


in  which 


(13e) 

(13f) 


F  = 


Similarly,  the  skew-symmetrical  part  of  the  problem  may  be  formulated 


by  having 

U,  V,  Oyy ,  Oyy,  0^2 

;  odd  in  z 

(14a) 

w’  Txz’  Tyz 

;  even  in  z 

(14b) 

to  be  continuous  across  the  plane  z  =  0  with  the  exception  of  the  crack 
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f 

t 

f 

{- 

region.  The  odd  quantities  in  eq.  (14a)  must  again  vanish  on  the  plane  of 
continuity.  On  the  crack  surfaces,  t  =  x  =  0.  Note  that  in  the  case 

At  Jr  fc 

of  skew-symmetry  azz  =  0  on  the  entire  plane  of  z  =  0,  i.e.. 


3  ( <j>  )  +  v  |  3<t>x  +  3(fy  = 

3Z  T^v  i  3x  3y  3Z 


(15) 


Eliminating  the  function  \p  by  means  of  eqs.  (7)  and  (15)  yields 


i  |  i*x  +  2  if* .  o . 

T-v  1  ax  3y  |  3z 


(16) 


Now,  let  g(x,y,z)  and  h(x,y,z)  be  two  harmonic  functions  related  to  <|>x  and 

♦y  ^ 


and  the  stresses  are 


%Q 

II 

2(l-v) 

A 

3y3Z 

2  36  4 

3Z 

■  z  s2g 
z  U? 

• 

(21a) 

II 

2(l-v) 

i!a_. 

3X3Z 

2—4 

3z 

■  z  s2g 

z  w 

• 

(21b) 

azz  _ 

2T“ 

32G 

z  izZ 

• 

(21c) 

ii 

-O-v) 

3_  ,ia 

3Z  v3y 

+  — ) 
3X' 

+  A. 

z  3X3y 

• 

( 21  d ) 

Txz  . 
2u 

(1-v) 

3  ,3£ 

W  'ay 

.  it)  4 
3X' 

3G 

•  M*  z 

32G 

3X3Z  * 

( 21  e) 

ii 

(1-v) 

_3 _  /  3h 

3X  '  3X 

■  SSL)  4 

3y 

■  3G.  z 
3y 

32G 

3y3Z  * 

( 21  f ) 

where 


(22) 


The  above  formulation  places  no  restriction  on  the  geometry  of  the 
planes  of  discontinuities.  However,  for  the  sake  of  definiteness,  the  sub¬ 
sequent  work  will  be  concerned  with  a  plane  crack  in  the  shape  of  an  ellipse. 

Temperature  Distribution  in  an  Infinite  Solid 

Consider  the  problem  of  an  infinite  solid  the  interior  of  which  is  ex¬ 
posed  to  uniform  temperatures  over  a  region  occupied  by  the  ellipse 

.  *-o  ■ 


Hence,  it  is  expedient  to  solve  this  problem  in  ellipsoidal  coordinates 
(C»n,c),  which  are  related  to  the  rectangular  coordinates  (x,y,z)  of  any 
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point  by  [9] 


a2(a2-b2)x2  =  (a^+s)  (a2+n)  (a2+c)  . 

(23a) 

b2(b2-a2)y2  =  (b2n)  (b2+n)  (b2n)  . 

(23b) 

a2b2z2  =  . 

(23c) 

where 

2  2 

<*>>£  »  0  >  n  ^  -b  >,  ;  -a  . 

In  the  plane  z  =  0,  the  inside  of  the  ellipse  is  given  by  £  =  0,  and  the  out¬ 
side  by  n  =  0. 

Let  the  temperature  distribution  T(x,y,z)  in  the  solid  to  be  an  even 
function  of  z  vanishing  at  infinity.  The  boundary  conditions  of  interest 
are 


T  ■  -T0  ,  5  =  0.  (24a) 

fl-O  ,  n  *  0  .  (24b) 


where  T  is  a  constant.  The  solution  of  this  problem  is  well  known  since 
o 

the  temperature  T(x,y,z)  is  equivalent  to  the  velocity  potential  of  a  per¬ 
fect  fluid  passing  through  an  elliptic  aperture  of  a  rigid  partition^*. 

Thus,  it  can  easily  be  shown  that 

T(x,y,z)  '  ‘  2KTkT  {  7$r =  'T°  W  •  (25) 


5See,  for  example,  [9],  p.  150,  eq.  (1). 
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Here,  K(k)  is  the  complete  elliptical  integral  of  the  first  kind  with 
argument  k2  *  1  -  (b/a)2  and 

Q(s)  =  s(a2+s)  (b2+s)  . 

The  variable  u  is  associated  with  the  Jacobian  elliptic  functions  snu, 
cnu,  — ,  and  should  not  be  confused  with  the  x-component  of  the  displace¬ 
ment  vector  u.  The  relationship  between  u  and  the  ellipsoidal  coordinate 
5  is 

5  =  a2  (2LU)2  =  a2  (sn"2u-l )  . 
sn  u 


From  eq.  (25),  the  temperature  gradient  may  be  computed: 


(26) 


On  the  plane  z  =  0,  eqs.  (25)  and  (26)  provide  the  correct  boundary  con¬ 
ditions  as 


Ho  »5  =  0 


T  = 


il 

9Z 


=  < 


Ui 


0  W 


=  0 


bfey 


(1  - 


x2  _  v2H/2 


a?  b 


9 


5  -  0 


0  ,  n  =  0  . 


Substituting  eq.  (25)  into  the  third  expression  in  eq.  (5)  gives 


9ij>  _  ,l+v\  aaTo  r  ds 

9i=  7tm 


(27) 


The  explicit  expression  of  ip(x,y,z)  is  not  required  for  this  part  of  the 
problem.  Terms  containing  the  derivatives  of  ^(x,y,z),  such  as  those  in 
eqs.  (12)  and  (13),  will  be  found  subsequently. 
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If  the  elliptical  region  of  exposure  is  maintained  at  a  uniform  tem¬ 
perature  gradient  Q0,  then  the  temperature  T(x,y,z)  is  an  odd  function  of 
z.  The  boundary  conditions  are 


|I=-Qo  .  5  =  0  . 


(28a) 


T  =  0  ,  n  =  0  . 


(28b) 


Aside  from  a  multiplying  constant,  T(x,y,z)  is  identical  with  the  velocity 
potential  for  axial  flow  past  an  elliptic  disk  in  an  infinite  fluid**.  From 
this  hydrodynamical  analogy,  it  is  found  that 


r"  ds 
2  s/Q(s) 


snu  dnu 
cn  u 


E(u)] 


.z  . 


(29) 


where 

u 

E(u)  =  /  dn^t  dt  . 
o 

and  E(k)  is  the  complete  elliptical  integral  of  the  second  kind.  Using  eqs. 
(28)  and  (5),  the  thermoelastic  potential  can  be  determined  from  the  con¬ 
ditions 


which  can  be  written  in  the  equivalent  form 


6See  [9],  p.  144,  eq.  (8). 
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♦  =  H(x,y)  ,4=0. 


-  q  n  =  o 

3  Z  * 


The  function  H(x,y)  may  be  considered  to  be  any  particular  solution  of  the 
equation 


32H  +  32H 
3x?  3y? 


fl+i aQn 
M-v'  2  * 


Furthermore,  if  H(x,y)  is  sufficiently  smooth,  a  suitable  solution  for  the 
thermoelastic  potential  may  be  taken  in  the  form  of  the  Newtonian  potential 
of  a  simple  layer  with  a  continuously  differentiable  density.  Without  going 
into  details,  it  can  be  verified  that  putting 


into  aT  =  2(j^-)  results  in  eq.  (29). 

Having  obtained  the  thermoelastic  potentials  given  by  eqs.  (27)  and 
(30),  a  more  detailed  treatment  of  the  thermoelastic  problem  is  in  order. 


The  Elliptical  Crack  Problem 

Let  the  surfaces  of  a  flat  elliptical  crack  in  an  infinite  solid  be 
opened  by  the  application  of  uniform  temperature  T0  and  temperature  gra¬ 
dient  Q0  as  described  in  the  previous  section.  In  the  absence  of  mechanical 
and  thermal  disturbances  at  sufficiently  large  distances  away  from  the 
crack,  the  displacements  and  stresses  are  assumed  to  vanish  as  z  approaches 
infinity.  From  the  knowledge  of  ^(x,y,z)  or  its  derivative  with  respect 
to  z,  the  complete  solution  of  the  present  problem  requires  the  evaluation 
of  one  function  f(x,y,z)  for  the  symmetric  part  and  two  functions  g(x,y,z) 
and  h(x,y,z)  for  the  skew-symmetric  part. 
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Suppose  that  the  prescribed  temperature  are  constant  across  the  upper 
surface  of  the  elliptical  crack  and  are  exactly  the  same  as  those  across 
the  lower  surface  given  by  eqs.  (24).  In  addition,  the  corresponding 
mechanical  conditions  are 

°iz  =  0,  £  =  0  ;  w  =  0,  n  =  0  . 

Thus,  eqs.  (12c)  and  (13c)  give 


obtained  from  eq.  (27).  Here, 


Upon  integration  with  respect  to  the  variable  z  yields 

94»  _  _  / 1 \  aaTn  f  X  dz 
9x  1-v  *  8K(k)  z  h2(a2+e)vt?(fT 
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(34) 


where  the  variable  of  integration  may  be  changed  from  z  to  5  by  the  re- 


Since  integrals  of  the  form 

00  00 

/  dz  /  x (s)  ds  . 
z  S 

may  be  evaluated  from  the  formula 

00  00  CD 

-/  dz  /  x(s)  ds  *  J  [z  -  z(s)]  x(s)  ds  . 
2  £  £ 


(35) 


(36) 


(37) 


(38) 
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eq.  (36)  takes  the  form 


/  If dz  =  O  •  ffffy  •*  Cz  1,-y 


rl+v.  aaT 


(39a) 


where 

I 


1 


I„  = 


-1/2 


ds 


•  1/2  ,  -3/2  «  -1/2 

/  s  (a2+s)  (b2+s)  ds  . 


In  a  similar  manner,  it  can  be  shown  that 


4lfdz  =  (pf>  -v  • 


(39b) 


where 


I3  ■  /"  (a2+s)'1/2  (b2+s)'3/2  t  1  -  4r  1  1/2  ds 

*  r  ac+s  b£+s 


I.  -  f  s1/2  (a2*sf 1/2  (b2+s)-3/2  ds  . 


The  displacements  in  the  plane  z  =  0  are  found  as 


8 


7The  evalu  . ion  of  eq.  (39a)  is  intimately  related  to  the  elliptical  punch 
problem  which  was  solved  only  in  part  by  Green  and  Sneddon  [8].  Their  solu¬ 
tion  can  be  completed  by  following  the  steps  going  from  eq.  (27)  to  (39). 
Incidentally,  a  complete  solution  to  the  elliptical  punch  problem  has  been 
given  by  Gal  in  [11]. 

®See  Appendix  (1)  for  the  results  of  integrals. 
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CE(u)  -  (|)2u  -  k2  . 


+  |^}[u-E(u)  +  k2  (40b) 


w  =  2(1 -v)  A  Lim  [  z.  /  -4L=  ] 
z-bO  C  Si/Q(s) 

=  iU^]A  Lim  [  z  fSniLdnu  .  E(u)]  |  .  (40c) 

ab^  z^O  ‘  cn  u  * 


Now,  let  uQ,  vQ,  and  wQ  be  the  displacements  of  the  crack  surfaces.  Then 
as  c  -*■  0: 

sn  u  -*  1 ,  cn  u  -►  0,  dn  u  -*•  j,  E(u)  ->  E(k)  . 


and  eqs.  (40)  simply  reduce  to 

Uo-.^EW  ,x  {  ,l-2v)  [Mil-1]  +  (1)2  m-1)  . 

'y  I  C(F»2  -  f^]  + 1  *  fjf}  )  • 

.  iLLvlA  0 . 4-  4)  1/2  . 

W0  ab  11  ?  b?' 


(41a) 

(41b) 

(41c) 


In  the  special  case  of  a  penny-shaped  crack,  a  =  b,  E(k)  =  K(k)  =  tt/2, 
eqs.  (41)  may  be  further  simplified  as 
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1  =  _  lit^SLk  £*J  t  W()  ,  2(H-v)acT0  (42) 

2  2  2 

where  p  =  r/a  and  r  =  x  +  y  .  The  normal  component  of  the  displacement 
w  of  the  crack  surface  aqree  with  eq.  (90)  in  [3]  only  when  the  integrals 
in  [3]  are  evaluated  properly**.  As  a  consequence,  the  variation  of  w0  with 
p  shown  by  Fig.  3  in  L3]  should  also  be  changed  accordingly. 

The  problem  of  a  uniform  steady  heat  flow  disturbed  by  an  insulated 
elliptical  crack  is  equivalent  to  the  one  of  assigning  uniform  temperature 
gradients  to  the  crack  surfaces.  The  thermal  conditions  are  shown  in  eqs. 

(28)  and 

Txz  =  Tyz  =  °’  £  '  0  :  U  '  V  '  °*  n  =  0  • 

must  be  satisfied  on  the  plane  z  =  0.  From  eqs.  (20a),  (20b),  (21e),  and 
(21f),  these  conditions  may  be  put  into  the  forms 

§+*!y  <§-£>  ’£  =  0-  <43a> 

*{  =  0-  (43b> 


9The  corrected  expressions  in  Olesiak  and  Sneddon's  paper  [3],  p.  253,  should 
read  as 

^(n)  =  2(1+v)a9oa  (Sin  n  ■  n  cos  n)  . 
it 


(88) 


and 

w(p)  *  [  0-p2)1/2  -  log  (1  +  ]  . 

Hence , 

w0  =  L0  (1-P2)1/2  .  (90) 

where 


and  e0  corresponds  to  T0  in  the  present  work. 

where  the  right  hand  side  of  eqs .  (43)  are  known  quantities  from  eq.  (30), 
and 

H*0,||-0,n»0  .  (44) 

Both  g(x,y,z)  and  h(x,y,z)  are  even  functions  of  z.  The  character  of  eqs. 

(43)  suggests  that  the  problem  can  be  readily  solved  by  putting10 

.  (45a) 

3h  ,  3_  3  .  . 

3Z  ^  3Z  3y^  0  •  (45b) 

Equations  (45a)  and  (45b)  represent  the  velocity  potential  of  a  rigid  el¬ 
liptical  disk  rotating  about  the  y-axis  and  x-axis,  respectively11. 

10If  H(x,y,z)  is  a  harmonic  function,  i.e.,  v2H(x,y,z)  =  0,  then  (x  ~ ■  - 

z  1_)  H  and  (y  L*  z  — )  H  are  also  harmonic. 

3x  Bz  By 

USee  [10],  p.  145. 
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By  virtue  of  eqs.  (45),  eqs.  (20)  render 


u  ■  -2(l-v)  X  ||  +  z  |p[  r  +  2(1-v)a  ]  , 
v  =  -2(l-v)  y  ||  +  z  |p-  [  r  +  2(1-v)n  ]  , 


(46a) 

(46b) 


w  =  2(1 -v)  *  -  (l-2v)  r  +  z  ||  . 


(46c) 


Equations  (21)  may  also  be  expressed  in  terms  of  the  functions  A(x,y,z)  and 
n(x,y,z) : 


_xx_ 

2y 


L.  ,  i_ 

az  vy  ay 


ar 


°zz 

2y 

Txy 

2y 


a2r 
z  az? 


a_ 

az 


3A 

3y 


an. 

axJ 


1*1  =  (l-v)  L(xli-yiE)+L(zL+l)r 
2y  ay  ay  ax  ax  az 

n 


■  az  + 

Z  [■ 

•  1L  + 

Z  [. 

az 

3*  . 

z  axay  *-r 

(z 

a_ , 

a2r 

ax? 


a2n_ 

3y?- 


a3 

+  (l-v) 


^3  °° 

a^7  [  /  dz  ] 


2y 


(l-v)  f-  (y  ^  *  x  |A)  +  f-  (Z  +  l)  r 

3X  3X  3y  3y  3z 


’  (1'V)  ^  C  /  Z(A'n)  dZ  ] 


(47a) 

(47b) 

(47c) 

(47d) 


(47e) 


(47f) 


where  the  following  contractions  has  been  made: 

aA  a n  an  92  00 

r  =  i|)  +  A  +  x9x  +  y  ay+z9z+  9X2  [  J  z(A-n)  dz  ] 


From  eqs.  (30),  (43),  and  (45)  follow  then 
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!?  <x !?  - 2  *  +  *!!7{xf  - ^r> +  [  0(A'n)  dz]l 


b\2 


3  .  3 

3Z  ^  3Z 


z  |y)  a +  v  (It  (y  If  - x  17)  - 13^7 1  Jz  z<fl-n> dz  ] ) 

-  <fe)  ffe  t 1  -  <?)2  §} ]  -y.«  =  °- 


Taking  A(x,y,z)  and  q (x,y,z)  in  the  forms 


[  ]  =  1/2  [  ]  (35?  +  B&I +  I-'  $ST  ’ 


the  constants  B  and  C  may  be  solved  from  eqs.  (48),  which  yield 

B  =  mlQ2  ~  m2Pl  c  =  mlP2  ~  m2Pl 
PlQ2  ’  P2<ll  *  dlP2  “  Q2P1 


The  constants  p^,  q^  and  riK  (j=l,2)  are 

P,  =  -  +  57^7  (  E(k)  -  K(k)  +  v  [(f)2  E(k)  -  K(k) 

-  (a2+b2^K(k)^-  2a2E(k)]  J  f 


q  -  ^  [  E(k)  -  K(k)  +  3  . 


m  = 


(]^)  .^[§}-u(|i: 


and 


P2  -  57^7  [  K(k)  -  (f)2  E(k)  +  (a2^b2^(k^-  2azE(k)  j  _ 


(48a) 


(48b) 


(49) 


(50) 
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q2  *  ^  (  k(k)  -  (i)2  E(k)  +  «  [  K(k)  -  E(k) 


+  (a2fb2)K(k)  -  2a2E(k)  ,  I 
a2  .  b2  1  I  ’ 


The  solution  is  essentially  complete  and  is  a  complicated  function  of  the 
material  constants  and  the  geometry  of  the  problem. 

A  straightforward  calculation  employing  eqs.  (46)  gives  the  displacement 
components1 ^  for  z  =  0 


Bx 


=  —2 ( 1 -v )  / 


v 

V.  J 


a3k2w 


Um  .  _L.  [  z  [ 
z-0  abz  ' 


snu  dnu 


cn  u 


-  E(u)  ]  ) 


(51a) 


Cy 


[  u-E(u)  ]  (|  (^)  .  (l-2v) 


K 

[  3B  +  (1  -  ^)  (B-C)  ]  ]  x  +  ^  l  . 

+  <  (j^)  -  (l-2v)  [  3B  +  2C  +  (^)2  (B-C)  ]j/(f)2  [  E(u)  -  k2  . 


snu  cnu 
dn  u 


u\ 


ja2k2  [  (lli)  -  (l-2v)  B  ]  +  (l-2v)  (B-C)  x2  J  u 

+  (l-2v)  (B-C)  x2  sn  u  cn  u  dn  u 


(51b) 


12The  results  of  the  integrals  are  tabulated  in  Appendix  (2). 
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In  the  same  way  as  before  by  setting  5  =  0,  the  displacements  of  the  crack 
surfaces  are  found: 


-  4(l-v) 
ab 


Bx 


Cy 


2  2 


1/2 


*  [  K(k)  -  E(k)  ] 


db^gQn 


(l-2v) 


(52a) 


[  3B  *  (1  -  y  (B-0  ]  ]  *2+  ]  • 

+  -  (l-2v)  [  3B  +  2C  +  (^.)2  (B-C)  ]j 

[  (|)?  E(k)  -  K(k)  ]  .y2 

-  (  a2k2  [  (]li)  -  (1-2v)  B  ]  +  (l-2u)  (B-C)  x2  )  K(k)  .  (52b) 


In  the  limit  as  a  =  b,  eqs.  (50)  shows  that 


B  ■  C  •  -  to  .  3#  • 

I -V  OTT 


Because  of  rotational  symmetry  of  the  circular  crack  problem,  the  displace- 

2  2  1/2 

ments  u  and  v  may  be  combined  as  (uQ  +  v£)  to  give  the  radia".  compo¬ 


nent  uf: 


3JF- 


(l+v)a^aQ0 


/l  -p^  ,  p  <  1 


(53) 
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13 

The  normal  displacement  w0  is 

•  (54a) 

and 

<&>  •  B&u  •  ST  (  ni-v)  p2  -  f  (2-V)  ]  sin-1  (1) 

-  (1-v)  I  p  »  1  .  (54b) 

Eqs.  (53)  and  (54)  check  with  eqs.  (26)  and  (27)  in  [5],  respectively,  when 
the  expression  for  w0  in  [5]  is  corrected.  The  temperature  gradient  Q0 
corresponds  to  t  in  [5]. 

The  computation  of  thermal  stresses  will  be  considered  in  the  next 
section. 


I^For  p  >  1 ,  w0  in  Florence  and 


Goodier's  paper  [5]  should  be  corrected. 


wp 

a2otQo 


(1 1.  i .  ]_) 

'2*2  2*p2' 


+  2(1-2v) 

45  nr 


(3  3 . 7 . j_\ 
22  Yp2 


where  2^1  (a,b;c;z)  is  Gauss'  hypergeometric  series.  Moreover,  these  hyper¬ 
geometric  series  can  be  reduced  to  elementary  functions  by  means  of  Gauss' 
recursion  formulas  [9]  and  some  properties  of  the  hypergeometric  series. 

The  results  are 


(1  3 .7.1_ 

2’2’2’p2 


(-)  -  /pz-1  ] 


There  is  also  a  misprint  in  [5],  eq.  (27),  where  the  factor  tt/4  should  read 
as  1/4.  .  26  - 


Crack-Border  Thermal  Stresses  and  Stress- Intensity  Factors 


A  knowledge  of  thermal  stresses  in  a  small  region  ahead  of  the  crack  is 
essential  in  the  investigation  of  the  stability  behavior  of  cracks.  For  in¬ 
stance,  the  Griffith-Irwin  theory  of  fracture  [7]  is  based  upon  the  concept 
that  the  onset  of  rapid  crack  extension  occurs  when  the  magnitude  of  crack- 
border  stress  field  or  simply  the  stress-intensity  factor  reaches  some  crit¬ 
ical  value.  For  a  given  configuration  of  the  crack,  this  value  will  in 
general,  depend  upon  the  properties  of  the  material  such  as  the  shear  modulus, 
Poisson's  ratio,  and  the  coefficient  of  thermal  expansion,  etc.  Hence,  at¬ 
tention  will  be  focused  on  the  determination  of  stress- intensity  factors 
kj(j=l,2,3),  where  each  one  of  the  k-factors  refers  to  a  particular  mode  of 
crack  surface  displacement.  The  elliptical  crack  problem  will  be  used  as  an 
example. 

The  "opening  mode"  of  crack  extension,  governed  by  k] ,  arises  in  the 
symmetrical  problem.  For  the  mere  purpose  of  finding  k-j ,  it  suffices  to  cal¬ 
culate  azz  in  the  plane  z  =  0,  outside  the  elliptical  crack  region.  By  way 
of  eqs.  (32)  and  (13c)14 


2  u  "  A^C1/2  (a2+c)VZ  (b2+C) 


tt? '  I 


"  (2s+a2+b2)ds 

5  ( a2+s ) ( b2+s ) /Q ( s ) 


_  00 

j_  /l+v\  aaTn  r  ds 

mV  •  4K(k)  \  ^rrr  * 


l+v  oTq_  j  J _  _ 

"  (  )  I  L  1/2 


ab‘ 


!-«■  2  '  E(k)  "  (a2+e)1/2  (b2+c) 


%  snu  cnu  _ 

17?  •  E(u)  +— -  ] 

dn  u 


+  u 
K[kT 


=  0  . 


(55) 


14Refer  to  Appendix  (1)  for  solutions  to  integrals. 
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In  the  neighborhood  of  the  crack  border,  the  ellipsoidal  coordinates  5 
and  c  have  the  limits^ 

-1/2 

S  =  2  abr  (a2  sin2*  +  b2  cos2*)  , 

c  =  -  (a2  sin2*  +  b2  cos2*) 


(56a) 

(56b) 


which  are  valid  for  n  =  0  and  r  «  1.  The  small  distance  r  in  the  plane 
z  =  0  is  measured  normal  to  the  crack  border  and  *  is  the  angle  in  the 
parametric  equations  of  an  ellipse.  Substituting  eq.  (56a)  into  (55),  the 
expansion  of  ozz  for  small  values  of  r  (or  s)  is 

°zz  =  (7^  F7 ^/^(a2  sin2<,)  +  1)2  cos **)  A“  +  0(r  )  •  (57) 

zz  1 -v  E(k )/  a  /2r 

The  stress -intensity  factor  k^  can  now  be  extracted  from  eq.  (57)  as 


The  variation  of  k^  with  the  angle  *  is  plotted  in  Fig.  1  for  v  =  1/3 
and  different  values  of  the  ratio  a/b.  A  glance  at  the  curves  shows  that 
k]  is  always  greatest  at  *  =  90°,  i.e.,  the  intersection  of  the  crack 
boundary  with  the  minor  axis  of  the  ellipse.  Thus,  crack  propagation,  if 
it  occurs,  would  first  take  place  at  the  point  (0,b,0),  and  tend  to  produce 
a  penny-shaped  crack.  This  is  a  brief  physical  interpretation  of  the 
Griff ith- Irwin  concepts  of  fracture. 


^See  [6],  eqs.  (52a)  and  (55). 
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The  remaining  strength  of  a  solid  whose  continuity  is  interrupted  by 
a  circular  crack  may  be  estimated  from  eqs.  (57)  and  (58)  by  setting  a  =  b 
and  £  =  r^  -  a^.  Consequently^ 

°zz  _  /1+v*  oTo  /  2  „  1 

27"  (T7}  w  (p  _1)  P>1  * 

and  (59) 

k,  =  (to  Ji  . 

I  I  -V  7T 


The  stress -intensity  factors  ^  and  k3  correspond,  respectively,  to  the 
"edge-sliding"  and  "tearing"  modes  of  fracture.  Their  evaluation  calls 
for  the  expressions  txz  and  TyZ  as  z->0.  It  is  found  from  eqs.  (47e),  (47f) 
and  (49)  that^ 


1/2 

1*2  .  [  ]  [  (l±v)  .^.Qq  +  B  +  vC  ]  /  _ 

2„  a?  -  b?  1-v  4E(k)  \  U2+s)MJT 


-  B  [ 


(2s+a2+b2)ds 
MU  \  (a2+s)(b2+s)/Q(s) 


-  / 


+  v/ 


ds 


] 


5  (b2+s)^IiT 

-Tjrr^h  2'  [  8  (1  -  -§-)  -  c  (^)1/2  0  -  — )  ] 

V'1  U-t)(b2-a2)  b2+$  a2+s  a2+c  b2n 


+  v  (B  -  C)  / 


s  ds 


(a2+s)(b2+s)/Q(s) 


,  n  =  0 


(60a) 


^The  factor  of  a/3  should  be  taken  out  of  eq.  (92)  in  [3]. 
l^The  integrals  are  tabulated  in  Appendix  (2). 
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o  o  V2 

L  iphtti .b£+<i  ]  .  |  [  (11*)  +  c  +  vB  ]  r” - * - 

2  ii  b2  -  a2  1  L  1-V  4E(k)  (  (b2+s)^T 

-  c  [  J_  -  f  <2s+a2+b2)ds _ t  v  r" _ ds _ ] 

v'QU)  5  <a2+s)  (b2+s)^Q(s )  $  (a2+s)/Q(s) 


¥¥¥? . t  b  <¥i)1/2  0  -  4-)  -  c  (¥i)'/,:  0  -  4-)  ] 


1/2 


Cl7^  ( c-c) (a2-b2)  b2+s 


a2n 


a2+C 


b2+c 


-  v  (B  -  C)  J 


s  ds 


{a2+s)(b2+s)*^Q(s) 


=  0 


(60b) 


where  B  and  C  are  given  by  eq.  (50).  Near  the  crack  boundary,  eqs.  (56)  may 
be  applied  to  simplify  eqs.  (60): 


Lcsljl/7  I  [  (i)2 

ab  ^  b  \  b 


D  r  1  t  v  b  sin  <j> 

B  -  t  J  [-x - 7-x - 5“ 

a^  sin^  +  b^  cos^<j> 


) 


-  B 


1/4 

(a2  sin2d,  +  b2  cos2^) _ +  Q  (rl/2) 

I  J2r 


111  -  2  sin  d)  fb  I  r  #bx2  r  R  ,  v  a2  cos  j>  x 
2  y  ab  U  l  L  V  J  V  sin2*  +  b2  cos2/ 


-  C 


(a2 


sin2j> 


b2  cos2<t> ) 


1/4 


/2r 


+  0 


(r1/2) 


(61a) 


(61b) 


Before  kj ( j-2 ,3)  can  be  found,  the  following  transformations 


xnz  Txz  S  ^  iy2  sin  8  . 

(62a) 

Ttz  -  ”  Txz  ^'n  ®  +  Tyz  ® 

(62b) 

should  be  performed  to  obtain  the  normal  and  tangential  components  of  the 
shear  stress.  In  eqs.  (62),  6  is  the  angle  between  the  x-axis  and  the  line 


normal  to  the  crack  border  in  the  plane  z  =  0.  It  is  not  difficult  to  show 
that  for  a  given  value  of  $,  i.e.,  a  fixed  position  on  the  crack  boundary, 

8  is  known  either  from 

-1/2 

sin  8  =  a  sin  <t>  (a2  sin2$  +  b2  cos2^ 
or 

cos  8  =  b  cos  <p  (a2  sin2<ji  +  b2  cos24>) 

Putting  eqs.  (61)  into  (62),  Tpz  and  xt  may  be  written  as 

Tnz  =  +  0  (r1^2)  ,  Ttz  =  |sL  +  0  (r^2)  . 


-  (a2  sin2<j)  +  b2  cos2^)  (B  cos2<j>  +  C  sin2<j>)  .  (63a) 

/at) 

* ^ ~v ^  (a2  sin2<ji  +  b2  cos24> )  (—  B  -  —  C)  sin  <j>  cos  <j>  .  (63b) 

/ab  b  a 

It  is  worthwhile  to  mention  that  in  the  skew-symmetrical  problem  of  an 
elliptical  crack,  there  exists  a  combination  of  "edge-sliding"  and  "tearing" 
movements  of  the  crack  surfaces  as  both  ^  and  k3  occur  simultaneously. 

Hence,  the  shape  to  which  the  crack  would  grow  will  depend  upon  a  function 
of  k^,  k^  and  is  no  longer  a  priori  evident  as  in  the  case  of  the  opening 
mode.  For  v  =  1/3,  the  values  of  ^  and  k3  against  <j>  are  shown  in  curves 
by  Figs.  2  and  3.  Note  that  k2  varies  with  <j>  in  a  manner  similar  to  k-j 
as  in  Fig.  1.  Inspection  of  rig.  3  reveals  the  interesting  fact  that  the 
four  points  at  which  the  crack  border  intersects  the  major  and  minor  axes 
are  under  the  action  of  "edge-sliding"  type  of  displacement  only,  since  k3 
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where 


vanishes  at  those  places.  The  maximum  values  of  k3  shift  as  the  ellipticity 
is  changed.  This  is  clearly  illustrated  in  Fig.  3. 

In  the  degenerate  case  of  decreasing  ellipticity,  i.e.,  as  b-*-a,  eqs . 

(63)  reduce  to 


(64) 


and  the  tearing"  mode  disappears  completely.  This  result  is  in  agreement 
with  that  obtained  by  Florence  and  Goodier  [5]^®. 

It  should  be  emphasized  that  the  angular  distribution  of  the  three- 
dimensional  thermal  stresses  near  thq_  crack  border  is  found  to_be  the  same 
as  those  in  an  elastic  body  undergoing  deformation  at  constant  temperatures 
[6].  The  two-dimensional  case  was  discussed  by  Si h  in  [12].  In  retrospect, 
this  justifies  the  application  of  the  Griff ith-Irwin  theory  of  fracture  to 
cracks  owing  to  thermal  disturbances. 


The  sum  of  the  normal  stresses,  (qr  +  oj  or  (oxx  +  ovv)  in  [5],  on  the 
plane  of  the  penny-shaped  crack  is  incorrect.  It  can  beyseen  from  eqs. 
(47a),  (47b)  and  (49)  in  the  present  paper  that 

(JV*  *°xx  +  Oyy\=0  ^  §C( k?°  +  *Btc'  1  O  '  \z  "  b?) 

M1M  (B  if  +  C  -V2  >5  =  0  . 


and  hence,  for  a  =  b,  the  above  equation  becomes 


■1-vx  +  °yy^z=o  _  1 


(  1  -V  ^  ~  ' 

1  +v  4  ij  aQQa 


1  [  (1-v)  +  1  (1+v)  _L„  ]  ,  p  <  1 

7T  3  A  -  0* 


which  is  the  revised  form  of  eq.  (28)  in  [5]. 
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Remarks  on  a  Class  of  Boundary  Problems 

A  more  general  approach  to  the  problem  of  finding  the  steady-state 
temperature  field  in  ellipsoidal  coordinates  will  be  considered.  Once  the 
structure  of  T(x,y,z)  is  known,  a  number  of  thermoelastic  problems  of 
cracks  may  be  solved  by  the  method  described  earlier.  For  brevity  sake, 
discussion  will  be  restricted  to  the  symmetrical  problem  where 

T(x,y,z)  =  T(x,y,  -  z)  . 

The  formulation  for  the  skew-symmetrical  temperature  problem  follows  in  the 
same  way. 

Suppose  that  the  surfaces  of  an  elliptical  crack  is  thermally  disturbed 
with  the  following  boundary  conditions: 

T  =  T(x,y)  ,5=0  .  (65a) 

Jj  =  0  ,  n  =  0  .  (65b) 

which  reduces  to  eqs.  (24)  for  T(x,y,0)  =  constant.  A  well-known  solution 
of  eq.  (1)  in  potential  theory  is  [13] 

T(x,y,z)  =  //  SlxliXlI  dx'dy'  ,  (66) 

E 

where 

-!  ?  ?  9  -1/2 
r  =  [  (x-x'r  +  (y-y'r  +  z2  ] 

and  q(x ' ,  y',  0)  is  the  density  of  a  distribution  of  mass  in  the  space 
(x1,  y',  z').  The  integral  in  eq.  (66)  is  extended  over  the  region  occupied 
by  the  ellipse  x^/a^  +  y2/b^  =  1.  The  requirement  is  that  T(x,y,z)  must  be 
sufficiently  regular  at  infinity.  This  is  satisfied  as 
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-►  00 


-1 

T  +  RQ  /  /  q ( X ’  ,y ' )  dx'dy '  ,  when  RQ 

z 

where 

o  o  o  V2 
Rq  =  (x2  +  y2  +  z2) 

Now,  eq.  (66)  may  be  inserted  into  the  last  of  eq.  (5)  and  the  result  can  be 
integrated  to  admit  the  representation 

^(x,y,z)  =  |  ^  /  /  q(x'.y')  log  (R  +  z)  dx'dy'  .  (67) 

Apart  from  the  factor  i£  (Jiii),  eq.  (67)  is  Boussinesq's  three-dimensional 

2  1  -v 

logarithmic  potential  for  a  disk  in  the  shape  of  z ,  whose  mass  density  is 
q(x,y) .  While  ij>(x,y,z)  has  continuous  derivatives  vanishing  at  infinity, 
the  function  \j>(x,y,z)  itself  is  unbounded  as  R  ->  ®  ,  since 


^  [  log  (R  +  z)  ]  .  /  /  q(x,y)  dxdy  ,  when  R  •*  ~ 

z 

The  boundary  condition,  eq.  (65a),  can  be  satisfied  by  taking 


T(x.y)  =  '-!  I  -atej&MS' 

2 7r  z  [  (x-x')2  +  (y-y')2  ]  '2 


(63) 


Integral  of  this  type  is  used  in  the  contact  problems  of  the  theory  of 
elasticity,  in  particular,  the  problem  of  a  perfectly  rigid  elliptical 
punch  on  an  elastic  half-space  [11]. 

In  ellipsoidal  coordinates,  the  temperature  field  may  be  taken  as 

T(x,y,z)  =  /  $(u>)  =£==  .  (69) 

C  "THs) 

where 


u)(  s)  =  1 
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and  $(u))  is  any  twice  differentiable  function  in  the  interval  (0,1)  with 
finite  one-sided  derivatives  at  the  boundary  points  of  the  interval.  Note 
that  eqs.  (65)  are  satisfied  by 


£[  = 

3z 


o  1/2  » 

-  —  UnO  /  *'(u) 

ab  ^ 


1/2 

ds  .  2  $(0)  [  n;  (a2+$)(b2n)  ] _ 

s4(sT  ab  (£-n)(£-c) 


.  (70a) 


and 


T(x,y)  =  j'*d--f--f) 

0  bz+s 


ds 

•Ws) 


(70b) 


O 

It  should  be  pointed  out  that  if  $(oj)  is  a  polynomial  of  degree  n  in  x 
and  y2,  then  T(x,y)  will  be  another  polynomial  of  the  same  degree  in  x2 
and  y2. 

Alternatively,  eq.  (69)  may  also  be  used  to  satisfy  the  boundary 
conditions 


f  q(x,y) 


3T  _ 


3Z 


=  < 


0 


,  K  =  o  . 

,  n  =  0  . 


By  taking  q(x,y)  to  be  a  function  of  the  variable 


Z  =  1 


x2  £ 

a2  "  b2 


(71a) 

(71b) 


then  eq.  (71a)  gives 

*(u))  =  -  ^  /“  (72) 

2  TT  0  /to-Z 

In  particular,  if 

q(Z)  -  q0  Zn  .  (73) 
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eq.  (72)  becomes 


*(w) 


n  +  1/2 

[  4s)  ] 


The  temperature  distribution  may  also  be  computed  as 


T(x,y,z)  =  -  f  [  u(s)  ] 

2A  r  (n  +  i)  r 


n  +  1/2  ds 


■WST 


(74) 


For  n  =  p  -  1/2,  where  p  is  a  positive  integer,  the  integrand  in  eq.  (74) 
can  be  expanded  and  the  integration  may  be  carried  out  term  by  term^. 

A  few  cases  of  eq.  (74)  are 


(a)  n  =  -  1/2,  eq.  (74)  reduces  to  eq.  (25)  with 

n  -  T0 

%  bK$y  • 

(b)  The  integral  corresponding  to  n  =  0  has  yet  to  be  evaluated. 

(c)  The  temperature  T(x,y,z)  for  n  =  1/2  is 


19 

In  the  special  case  of  a  =  b,  eq.  (74)  can  be  integrated  for  all  values 
of  n  giving 


T=^^inil)  (k,) 

JiT  r  (n  +  |-) 


2(n  +  1 ) 


I 


2(n  +  1 )  »  S  =  0  • 


where 


K 


dt 


o  (dn  t) 


2n  * 
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»  e  =  o  . 


Conclusion 

Three-dimensional  problems  of  cracks  in  thermoelasticity  have  been  formu¬ 
lated.  The  displacements  and  stresses  are  expressed  in  terms  of  potential 
functions,  which  are  valid  for  planes  of  discontinuities  bounded  by  smooth 
curves.  Of  fundamental  interest  is  the  problem  of  a  flat  elliptical  crack 
whose  surfaces  are  thermally  disturbed.  Exact  solutions  are  obtained  and 
reduced  to  the  limiting  cases  discussed  by  previous  authors  [3-5]. 

Several  essential  features  of  the  stability  behavior  of  an  elliptical 
crack,  heated  on  its  surfaces,  are  discussed  in  connection  with  the  Griffith- 
Irwin  theory  of  fracture.  The  concept  of  stress-intensity  factor  is  intro¬ 
duced.  In  general,  all  three  types  of  k-factors  are  present  in  structural 
members  undergoing  thermal  changes.  Therefore,  a  criterion  of  fracture 
will  have  to  depend  on  a  function  of  k-| ,  k2,  and  kg  to  reach  some  critical 
value,  say  fcr  at  which  unstable  crack  extension  would  start.  The  function 
fcr  must  be  determined  experimentally. 

APPENDIX 


(1)  The  folloving  integrals  are  essential  in  the  calculation  of  dis¬ 
placements  and  thermal  stresses: 


f  ds  =  1 

J  /A  l  .  \  •  U  , 


c  ARsT  a 


(75a) 


/ 


ds _ .  2 


£  (a2+s)/Q(s)  a2k2 


/ 


ds  .2 


r  (b2+s)^TiT  a3k2  'b 


[  u  -  E(u)  ]  , 


[  (f)  E(u)  -  u  - 


^a2-b2^  sn  u  cn  u  j 


b2 


dn  u 


(75b) 


(75c) 
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00 


J  --**-£=  •  -ig  [  E(u)  -  k'2u  ]  , 

5  (a2+s)/Q(s)  ak2 


(75d) 


r“  s  ds  .  2  c  r/..i  ,.2  sn  u  cn  u  T 

J£  dAsl^lsT  '  ak7  [  E(U)  +  k  — dn  u  ]  ’ 


(75e) 


f _ ds _ .  _1 —  [  [iV|  E(u)  .  2  u  -  («1£)  sn  u  C",  a  ]  , 

£  (a2+s)(b2+s)»/t)(s)  a5k4  b2  b2  dn  u 

( 75f ) 


J  — [  (1  +  E(u)  -  2  u  -  sn  u  cn  u  . 
f  (a2+s)2/Q(s)  3  a3k2  k2  k2 

.  dn  u  ]  .  (75g) 


s  ds 


.  2 


2.  .2 


£  (a2+s)(b2+s)/Q(s)  a3k2  a2-b2 


[  <*=%>  U  ■  j-  E(u)  ♦  s".Acn_u,  J  .  (75h) 


dn  u 


(2)  For  the  skew-symmetrical  part  of  the  problem,  it  is  necessary  to 


find 


00  00 

/  z  .  A  (x,y,z)  dz  ,  /  z  .  n  (x,y,z)  dz 

z  z 

where  the  functions  A  (x,y,z)  and  «  (x,y,z)  are  themselves  integrals  given 
by  eg.  (49).  Hence,  integral  of  the  type 


/  ( 


a2+s  b2+s 


+  *--  1) 
s 


/Q(s) 


]  dz 


is  involved.  Using  formula  given  by  eq.  (38),  eq.  (76)  becomes 

i .  .1/"  (jsL  +  -iL  +  £.  D2 .  s_ds  _ 

8  £  a2+s  b2+s  s  ^(s) 


(77) 


From  eq.  (77)  follows  then 

il  =  .  1  f  ( JiL  +  J—  +  11-1)  _ _ x2  j°° _ sjs _ 

3x2  2  ^  a2+s  b2+s  s  (a2+s)v/Q(s)  £  (a2+s)2/Q(s) 

(78a) 
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00 


33I  _ 


3x23y 


=  -y  / 


s  ds 


2  xzy  g 


3/2 


33I  _ 


t  (a2+s)(b2+s)»^IsT  (u2n)3/2(b2H)1/2U-n)U-s) 

3/2 


3xay< 


=  -  x  / 


s  ds 


2  *y2  c 


5  (a2+s)(b2+s)^?T  ^nJ’^+eP^U-bJU-c) 


.  (78b) 


.  (78c) 


Similar  expressions  for  partial  derivatives  of  I  with  respect  to  x,y,z 
may  be  obtained  for  the  computation  of  stresses  and  displacements.  It  is 
now  apparent  that  the  solutions  to  the  integrals  in  eqs.  (78)  can  be  ob¬ 
tained  from  eqs.  (75). 
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